Astronomy & Astrophysics manuscript no. coher 


© ESO 2009 


August 28, 2009 





Q^; Solar active regions: a nonparametric statistical analysis 

O ; 

I J. Pelt^, M. J. Korpi^, and I. Tuominen^ 

^ , ^ Tartu Observatory, 61602 Toravere, Estonia 
<^ ■ ^ Observatory, PO Box 14, FI-00014 University of Helsinki, Finland 

Received — ; accepted — 

ABSTRACT 



00 



00 



I 

o 

5h 



Context. The sunspots and other solar activity indicators tend to cluster on the surface of the Sun. These clusters very often 
occur at certain longitudes that persist in time. It is of general interest to find new and simple ways to characterize the observed 
distributions of different indicators and their behaviour in time. 

Aims. In the present work we use Greenwich sunspot data to evaluate statistical but not totally coherent stability of sunspot 
distribution along latitudes as well as longitudes. The aim was to obtain information on the longitudinal distribution of the 
underlying spot-generating mechanism rather than on the distribution and migration of sunspots or sunspot groups on the solar 
C/3 , surface. Therefore only sunspot groups were included in the analysis, and only the time of their first appearance was used. 
, ^ , Methods. We use simple nonparametric approach to reveal sunspot migration patterns and their persistency. 

Results. Our analysis shows that regions where spots are generated tend to rotate differentially as the spots and spot groups 
^—i themselves do. The activity areas, however, tend to break down relatively fast, during 7-15 solar rotations. 
^ Conclusions. This study provides a challenge for solar dynamo models, as our results are consistent with the presence of a non- 
, axisymmetric spot-generating mechanism experiencing differential rotation (known as phase mixing in dynamo theory). The 
■ new nonparametric method introduced here, completely independent of the choice of the longitudinal distribution of sunspots, 
' was found to be a very powerful tool for spatio-temporal analysis of surface features. 
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' barrassment of riches" . When spatio-temporal proper 

■ ■ ties of the smaher features - say spots, flares etc - are The answers to these kind of questions depend very 

treated with weU-estabhshed vigour, the analysis of spa- much on the method of analysis used. Very often subjec- 

tially larger or temporarily longer patterns is very com- tive judgement is involved, either through steps of visual 

plicated. Even the nomenclature of the phenomena is processing or through involvement of freely chosen proce- 

not weU-estabhshed - for instance the time-space clus- dure parameters (bin sizes, zone widths, detection limits, 

ter of the local phenomena can be called as "active preselection criteria etc). 

longitu de" (Lo sh[19381 Vitinskij [MS]) , " Sonne flecke nherd " Yrom the statistical analysis point of view we can di- 

(Becker |1955p , "active region" (Bumba & Howard |1965p , vide the previously used methods along two lines: how the 

"sunspot nest" (Castenm iller et al.|1986p , "comple x of ac - input data is transformed before computing final statistics 

tivity" (Gaizauskas et al. |1983p or "hot spot" (Bai |1988p . and what kind of statistics are used. Some typical but ran- 

There is, in addition, a problem with proper definition of dom examples: 
such extended patterns. 

It is generally thought that the tracers of solar activity 

- sunpots, fiares etc - are randomly generated manifesta- — Aggregated data (daily Wolf numbers) and correlation 

tions of the larger scale mean magnetic field of the Sun analysis (Bogart [1982p . 

generated by a hydromagnetic dynamo process. An anal- — Raw heliographic longitudes and longitude- wise bin- 

ogy with a submerged animal blowing out bubbles is quite ning (Trotter & Billings [19621 Warwick I1965p . 

— Transformed (using trial rotation velocity) longitudes 

Send offprint requests to: J. Pelt and statistic (Bai ll987| . 



1. Introduction appropriate in this context (see Bai 12003"]) . What can we 

tell about the swimming speed and size of the animal, if 
Modern observations of the Sun are so rich in de- ^^^^ ^^^^^^ ^^^^y^^ observable? How deep in water 
tail that astronomers are eventually stroken by "em- ^j^^ animal? 
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— Transformed (using latitude-dependent rotation veloc- 
ities) longitudes and pattern matching (Usoskin et 
al. duUEi Pelt et alHOOi hereafter PBKT), 

— Spherical harmonic decomposition and time series 
analysis of mode amplitudes, phases and phase-walks 
(JuckettEOOa]). 

The third important aspect of the analysis is time cover- 
age of the observations. It is quite easy to find recurrent 
patterns in short time series, but coherence tends to break 
down very fast for longer datasets. 

In this paper we in a certain sense try to return to 
the square one, back to the very basics. Using very simple 
considerations and avoiding all freely chosen parameters 
we try to get answers to the questions: 

— Is there a tendency for surface elements to occur at 
certain longitudes that persist over time? 

— How this persistence and differential rotation of the 
surface elements are connected? 

— How long typical correlations in activity persist? 

Our aim is not so much to perform another statistical 
analysis of the well-known and already extensively anal- 
ysed data, but to introduce a new nonparametric method 
of analysis involving no physical, geometrical or statistical 
prior assumptions. In Section [2] we introduce our method 
of analysis, in Section[3]we present the results obtained for 
the Greenwich sunspot data set, and finally in Sections 2] 
and [5] we discuss our results in the light of previous sta- 
tistical analyses. 

2. Method of analysis 



2.1. Nonparametric method 



.(1) 



Let us assume that we have two sets of longitudes: X\ ' ,i 

(2) 

1, . . . , iV and ,.j — 1, • • ■ , M. Their values belong to 
the interval 0° < A < 360° and we assume that N <= M 
(if otherwise, we can always swap the sets). We want to 
characterize the similarity or the difference between the 
two longitude distributions somehow. The general theory 
of directional measurements is considered in mathemat- 
ical statistics (see for instance the latest monograph by 
Mardia & Jupp 120001 and references therein) , but here we 
need a more specific method, namely one without any un- 
derlying statistical assumptions or parametric models for 
the distributions involved. 

We propose the following very simple nonparametric 
method: The circular distance between two longitudes 
and A; we define as usually done 

AAfc,i =min(|Afc-An,360°-|Afe-A,|). (1) 

Let us take a particular longitude A^^"* from the first set. 
Among the longitudes of the second set there is always a 
value whose circular distance from the selected value is the 
smallest, let us denote this distance as AAfc. All together 
we can compute N such values - for each longitude in the 



first set. Now we select the longest distance among them 
and denote it simply as A. It is quite clear that, in the 
particular case when the first set is just a subset of the 
second one, A = 0. If the sets differ, then A > 0. In 
principle such max-min distance between two longitude 
sets is already a useful statistics; its full power, however, 
is revealed if we properly normalize it. 

For a particular set sizes N and M we can compute 
the mathematical expectation of A for completely random 
distributions of longitudes in both sets. Let us denote this 
expectation as A. Our final statistics which measures sta- 
tistical distance between the two sets of longitudes is then: 



-I- 



(2) 



If we want to stress that the distance D is computed for 
two particular indexed longitude sets, say for index n {N 
longitudes) and index m (M longitudes) then we use no- 
tation 



D{n, m) 



A. 



A 



N,M 



(3) 



The mathematical expectations A.n.m depend only on the 
integers N and Af, and can be pretabulated. In our calcu- 
lations we used approximations obtained from randomly 
generated longitudes for 10000 statistically independent 
runs. 

It is quite obvious that for absolutely random pairs of 
longitude sets our distance will have a value around 1 . For 
weakly correlated sets values are less than one, and values 
higher than 1 can occur when the longitude sets involved 
are constrained in a certain way due to which they cannot 
form all the patterns which occur for randomly generated 
sets. In the case of sunspot groups, for instance, the dis- 
tributions are constrained by group sizes. Randomly gen- 
erated points can fall arbitrarily close, which is not true 
for sunspot longitudes, because for them the group centres 
are separated by definition. 

Having now the statistic to measure distances between 
different distributions of longitudes we can go further. For 
a sequence of longitude sets we can compute a mean dis- 
tance between neighbouring sets: 



D 



k=l 

K - 1 



(4) 



where K is the number of the sets. We can also investi- 
gate how the distance depends on the mutual positions of 
particular sets: 



C{1) 



'j^ D{k,k + l) 

k=l 

K -I 



(5) 



Eventually D = (7(1). The statistic D allows us to in- 
vestigate rotational properties of the sunspot groups and 
statistic C{1) will be used to estimate how persistent the 
longitudinal correlations are. 
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2.2. Rotation and frames 

Heliographic longitudes are defined using the so called 
Carrington frame, which rotates against fixed stars with 
the exact period of Pc = 25.38 days. The mean rotation 
period if observed from the Earth is Pq = 27.2753 days. 
The Carrington frame is a formal construct and real fea- 
tures on the Sun need not to follow it exactly. 

Let us fix a certain longitude A*^*^-* of a particular per- 
sistent feature on the Sun rotating with the Carrington an- 
gular velocity. Then its longitude for different Carrington 
rotations i will be fixed: A,[^^ = [A^^^ + i x 360°] = A^'^) , 
angular brackets denoting here and below reduction to 
the interval (0,360°). Because the angular velocity of the 
Carrington frame is flc = degrees per day we can 
rewrite cycle dependent sequence of longitudes as 



= 0,1,, 



(6) 



The actual angular velocity of an arbitrary feature on the 
Sun need not to be exactly ilc- Let longitude of the first 
occurrence of such feature be A. Then cyclic reoccurrences 
of it can be described using a correcting term AQ: 



A,; 



[x + ix (nc + An)Pc],i^o,i, 



(7) 



The corrected frame rotates against Carrington frame 
with angular velocity Afi degrees per day. For convenience 
we introduce also notion for siderial angular velocity of the 
accelerated or decelerated frames = fie + AfJ. 

In what follows, we measure the angular velocity in 
degrees per day, latitude in degrees, and periods in days, 
and give the values in these units. 

2.3. Algorithms 

All ingredients of the method of analysis described, we can 
now formulate our basic algorithms. 

As an input data we use a set of time tagged longi- 
tudes ti,Xi,l — 1,...,L, amounting to L pairs of data. 
Using the time points ti we divide records into subinter- 
vals with the length 27.2753 (Carrington rotations). This 
procedure is not absolutely exact because the observation 
timing depends on the somewhat excentric orbit of the 
Earth. Fortunately the errors involved are small and we 
can ignore them. From the point of generality and ob- 
jectivity our choice is quite natural. Historical observa- 
tions are all done from Earth and consequently the fea- 
tures can be observed only half a time. However, during 
the Carrington rotation we can record what happens at 
all longitudes. As far as timing is considered, due to the 
rotation some processes can actually start earlier than ob- 
served. This excludes short-living processes (shorter than 
Carrington rotation) from our analysis. 

It is also possible to divide observations into longer 
subintervals. Then we increase statistical stability of our 
estimates (more observations in subsets) but loose resolu- 
tion in time. We consider time step with the length of one 
Carrington rotation to be optimal. 



We assume that the features on the surface of the Sun 
rotate with angular velocity which is different from the 
Carrington velocity Qc- For a certain trial angular veloc- 
ity Q and for each Carrington cycle i we can compute 
longitude corrections: 



A, 



i X AnPc ^ixifl- nc)Pc,i = 0, 1, . 



(8) 



By substracting rotation number dependent corrections 
from measured longitudes and properly reducing results 
to interval (0, 360) we build transformed longitudes: 



A(^) = [A,-A,] 



(9) 



They can be analysed using the statistics introduced 
above. We can also say that we transform longitudes in the 
Carrington frame into longitudes in the comoving frame. 
The frame rotation velocity is a free parameter of the 
procedure. We expect that the distributions of the trans- 
formed longitudes depend on 11 and the highest level of 
correlation in the longitude distribution will show up as a 
minimum of the distance statistic D. 

First we compute how the mean distance between 
neighbouring rotations D depends on angular velocity fi. 
Then we can use the best value (producing the highest 
level of correlation) for angular velocity to compute how 
distances depend on the interval between rotations (using 
the statistic C{1)). 

3. Data analysis 

Here we describe how we apply the presented statistical 
method to study the particular case of sunspots. 

The most comprehensive (in time) compilation of 
sunspot data was downloaded from the Science at NASA 
web sit£0. The same minor corrections as in PBKT were 
introduced. In this paper we used the full data set cover- 
ing years 1874-2008, or in terms of Carrington rotations, 
the rotations 275-2074. From all the data base records we 
chose only sunspot groups, leaving out single spots. In 
this way all the entries in the final set have equal statis- 
tical weight. For each sunspot group we selected only the 
record of its first occurence. This is important aspect of 
our analysis. We do not track sunspots as they rotate, but 
are interested in the movement of the underlying spot- 
generating structures. The final compiled data sets cover 
rotations 275-2074 with 16053 records for the Northern 
hemisphere and 275-2071 rotations with 15858 records for 
the Southern hemisphere; the compiled data sets are avail- 
able on the we 



3.1. Mean angular velocity 

For the first approximation we can assume that the mech- 
anism generating the sunspots rotates as a rigid body. 
Then we can measure its angular velocity using D statis- 
tic by comparing different longitude correction schemes. 



^ http:/ /solarscience. msfc.nasa.gov/greenwch.shtml 
http:/ /www.aai.ee/~pelt/soft.htm 
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n (deg/day) a (deg/day) 



Fig. 1. Statistic D{n) for the full dataset of the Northern 
hemisphere (thick line) and for the full dataset of the 
Southern hemisphere (thin line) . Minima indicate the best 
fitting comoving velocities. 



and choosing the one that produces the lowest value of 
the statistic. The results can be best illustrated by dis- 
playing D as a function of il - the actual siderial angular 
velocity of the frame. In Fig.[T]such functions are displayed 
for both solar hemispheres. As we see both curves show 
very clear and indicative minima. The absolute minini- 
mum for the Northern data is positioned at 14.348 and 
for Southern data at 14.403. The curves themselves are 
somewhat fluctuating so that we found useful to estimate 
the minima using local fits of the fifth-degree polynomials 
also. The resulting values from the fitting procedure are 
14.365 ± 0.001 and 14.409 ± 0.002 for North and South, 
deviating from the absolute minimum values by 0.017 and 
0.006, respectively. 

From these results we can see that the Southern part of 
the mechanism tends to rotate slightly faster. However the 
difference is too small, especially if to take into account 
the roughness of our method, to be conclusive. 

3.2. Differential rotation 

Sunspots and other activity indicators rotate with differ- 
ent angular velocities at different latitudes. By tracking 
particular objects in time it is possible to build a smooth 
curve to reveal the overall pattern of such differential, lat- 
itude dependent, rotation. Our statistic Z)(il) does not 
track single sunspots or the actual movement of sunspot 
groups, as we include only the first appearance of the 
sunspot groups. This way we can check whether the spot- 
generating mechanism itself rotates differentially or not. 
For that purpose we divided the observed groups into four 
subsets along latitudes (per hemisphere) and computed 
D{^) for every group. The latitude limits for the subsets 
where chosen to make them as equal in size as possible. 
The typical curves are shown in Fig. [51 The exact deter- 
mination of the minima for the curves is somewhat com- 
plicated. If we locally fit polynomials into the curves as 
we did above, we can get estimates with high formal pre- 



Fig. 2. Statistic D{n) for the latitude strip 0.0 - 9.8. 
Thick line - Northern hemisphere, thin line Southern 
hemisphere. Minima indicate the best fitting comoving ve- 
locities. 
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Latitude 
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Angular 


7b d 




range 




velocity 




North 


0.0-9.8 


4013 


14.430 


80.1 




9.9-14.4 


4015 


14.224 


71.8 




14.5-20.0 


4048 


14.371 


42.4 




20.1-60.0 


3977 


13.974 


3.7 


South 


0.0-9.8 


4095 


14.496 


78.3 




9.9-14.4 


3847 


14.423 


74.8 




14.5-20.0 


4143 


14.370 


41.3 




20.1-60.0 


3809 


13.964 


3.9 



Table 1. Differential rotation 
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Fig. 3. Differential rotation curve. Thick curve - fit for 
Southern data, thin curve - full fit. 



cision (0.001-0.002). The differences between the absolute 
numerical minima and the fitted minima, however, can 
be quite large (up to 0.045). Therefore we can claim that 
the probable statistical errors of the obtained minima are 
around 0.02 degrees per day. 

The full set of the absolute minima for all the eight 
curves is given in Table [221 In Fig. [31 the traditional least 
squares fit of the obtained angular velocities f2i is pre- 
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sented. We estimated parameters A and B for a simple 
model: 

n,=A + Bsm'^{0,), (10) 

where ^^-s are mid-latitudes of the eight belts. The es- 
timated parameters were A = 14.438 ± 0.035 and B = 
— 1.13 ± 0.17 correspondingly. For the Southern hemi- 
sphere the latitude dependence of the angular velocity is 
monotonically decreasing polewards and resembles quite 
well the curves obtained from sunspot tracking. For the 
Northern hemisphere the latitudinal behaviour is not 
monotonous, and the two latitude bands nearest to the 
equator are somewhat slower than expected. If we used 
only Southern velocities for the fitting procedure, the re- 
sulting fit was better with parameters A — 14.490 ± 0.001 
and B = —1.276 ± 0.044. At the moment we do not try to 
find any statistical or physical interpretation of this result; 
more importantly to us, the results obtained in this sec- 
tion clearly demonstrate that the simple nonparametric 
method can be succesfully used to study the differential 
rotation of the solar activity tracers. 

3.3. Break down times 

The results of the previous section clearly show that 
the longitudinally concentrated spot-generating mecha- 
nism is subject to differential rotation. Kinematic mean- 
field dynamo theory predicts (e.g. Krause&Radler 1980) 
that in the parameter regime where nonaxisymmetric dy- 
namo modes can be excited, the nonaxisymmetric modes 
are non-oscillatory and rotate rigidly with angular veloc- 
ity different from the overall rotation period. The phe- 
nomenon of phase- mixing, i.e. the nonaxisymmetric modes 
becoming affected by differential rotation, is against these 
predictions; our results, however, are consistent with the 
phase-mixing effect. 

Let us now try to quantify the effect of the differential 
rotation on the nonaxisymmetric structures by calculating 
the characteristic time needed to break down a longitudi- 
nally elongated stucture for different latitude strips. Using 
the estimated B values from Eq. [TO]we can define a break- 
down time for the strip of latitudes (6*1, 62): 

B(sin Oi — sm B2) 

where A0 is the phase distance over which the hypotheti- 
cal longitudinal pattern can be regarded to be destroyed. 
A reasonable value for the parameter A(/) comes from a 
following simple observation. Let us assume that at the 
latitudes Oi and O2 we have K observations. Let the longi- 
tudes coincide for the starting point in time. The statistics 
introduced above are based on finding the nearest "neigh- 
bours" . In our case for every observation at one latidude 
there is exactly one "neighbour" at the other. To break 
the ties between the neighbours we need to rotate the ob- 
servation at the other latitude. That means we need to 
have relative shifts which are longer than half the dis- 
tance between two consequtive observations or formally 
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Fig. 4. Statistic C{1) computed for the Northern latitude 
strip 0.0 — 9.8. Horizontal line is the mean value of the 
last 25 points. It indicates approximate asymptotic level 
of correlation. 
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Fig. 5. Statistic C{1) computed for the Northern latitude 
strip 9.9 - 14.4. 

A(/) = 0.5 X 360°/^, where is a mean number of obser- 
vations per Carrington rotation. Applying all this to the 
same latitude intervals as in Sect. 13.21 we obtain the break 
up times listed in Table [321 

3.4. Decorrelation time 

So far we have demonstrated that the sunspot group dis- 
tributions along longitudes for sequential Carrington ro- 
tations are correlated. We also computed the best fitting 
mean angular velocity Q for several latitude strips, and 
as a result found out a clear differential rotation pattern. 
Next we are interested in estimating the approximate life- 
times of the correlated features found from the sunspot 
data. For that purpose we use the obtained mean angular 
velocities for each latitude strip and compute C{1) curves 
to learn how fast the correlation between appropriately 
rotated longitude sets fades off. 

The results of this kind of analysis are displayed on 
Figs. |4]to[TlJ On these plots the thick lines indicate mean 
level over the last 25 points, which gives an asymptotic 
level of convergence. As can well be seen from the plots. 
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Fig. 6. Statistic C{1) computed for the Northern latitude Fig. 9. Statistic C{1) computed for the Southern latitude 
strip 14.5 - 20.0. strip 9.9 - 14.4. 
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Fig. 7. Statistic C{1) computed for the Northern latitude Fig. 10. Statistic C{1) computed for the Southern latitude 
strip 20.1 - 60.0. strip 14.5 - 20.0. 
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Fig. 8. Statistic C{1) computed for the Southern latitude Fig. 11. Statistic C{1) computed for the Southern latitude 
strip 0.0 - 9.8. strip 20.1 - 60.0. 



the inherent scatter of the C[l) curves is quite large; this 
is due to the physical variability of the activity level and 
the roughness of the statistic. Therefore it is hard to fix 
the point where the asymptotic level is achieved. Some as- 
pects of the curves, however, are quite indicative. First - 
the shortest decorrelation time is obtained for the highest 
latitudes. This is obviously because the width (in degrees) 
of the high altitude strips is wider and covers very differ- 



ently rotating spot groups. Secondly - the strips nearest 
to the equator show the longest correlations. We can quite 
safely claim that certain level of correlation is visible up 
to time span of 15 rotations. 

Comparing the estimated decorrelation times for dif- 
ferent latitude strips with the break-down times from 
Table 13.21 we can see that for lower latitudes the decorre- 
lation times obtained from our analysis are much shorter 
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than the estimated break-down times. For the highest lat- 
itude strip the decorrelation time is of the same order of 
magnitude than the break-down time. Part of this effect 
could be due to the enhanced diffusion of the field due to 
the stretching by the shear in angular velocity; the latitude 
dependence (shorter decorrelation times at high latitudes 
with the largest relative shear) would support this inter- 
pretation. It seems likely, however, that at least at lower 
latitudes, stretching and enhanced turbulent diffusion act- 
ing on the magnetic field due to differential rotation are 
not the only effects at play. 

4. Discussion 

To put our results into general context we will compare 
them with a sample of previous analyses. 

Large amount of the solar variability research is not 
based on a full set of sunspot observations but on some 
aggregated form of data. Most typically the daily Wolf 
sunspot numbers are used. For instance Bogart (1982) 
analysed these numbers using autocorrelation functions 
and power spectra. The major results were quite similar to 
ours - the rotation period around 27 days was detected and 
persistence of activity zones was claimed to be of the or- 
der of 10 solar rotations. In principle correlation functions 
and power spectra can be considered to be parameter-free 
statistics. The aggregated nature of the Wolf numbers, 
however, does not allow to analyse latitude dependence of 
the active clusters. 

There is a number of analyses which use longitudinal 
phase binning of the surface features. For instance in a 
series of papers Bai (|1987I1988|) used comoving frames (as 
in our work) to seek rotation velocities which enhance sta- 
tistical contrast of longitudinal distribution of solar flares. 
The transformed longitudes for each trial rotation veloc- 
ity were binned into 12 bins and variance of obtained dis- 
tributions was computed. The possibility of differential 
rotation was not taken into account. To study the per- 
sistence of particular active regions were visually tracked 
and displayed as "family trees" . In describing his results 
the author proposed a general scheme to characterize hi- 
erarchical patterns of solar activity: 

— single events (sunpots, flares), 

— active regions, 

— activity complexes, 

— active zones. 

The lifetimes of the activity centres increase with hierar- 
chy - from days to several years. In a later work Bai pOOSp 
used Rayleigh-type statistics to analyse transformed lon- 
gitudes. He computed standard spectra which are sensitive 
to unimodal distributions and spectra which are sensitive 
to bimodal distributions. As a result he found that some 
characteristics of the longitude distributions are rather 
persistent in time, even up to decades. Our results describe 
average behaviour of the solar activity and consequently 
some long-lived elements do not have strong influence, as 
they are mixed with other elements whose lifetimes are 



shorter. It should be also stressed that D and C{1) statis- 
tics do not depend on any assumption about modality of 
the underlying variability. All the "modes" are automati- 
cally accounted for. 

Probably the most popular method to study the kine- 
matics of the solar surface features is a standard power 
spectrum analysis and its variants (just an example - 
Temmer et al. [2004| Giordano [2008]) . This kind of anal- 
ysis can be applied to latitude strips and in this way the 
differential rotation can be taken into account. From the 
first sight Fourier analysis seems to be essentially nonpara- 
metric. However, the fact that it uses single harmonics as 
base functions prescribes certain form of preferred activ- 
ity distributions. The results of the Fourier methods are 
often given as a list of certain periods which show up in 
power spectra or on wavelet plots. The periodicity claim 
itself is quite a strong statement, as it is often very hard to 
find physically solid timing mechanisms for periods which 
strongly differ from the obvious one - that of solar rota- 
tion. 

We want to stress here that in the proposed statistical 
method no assumptions about the particular form of the 
activity indicator distributions are made. Even more - the 
statistics D and C{1) are not seeking certain clusters or 
other kind of patterns, they are just used to check whether 
the "birth places" of surface elements are correlated or 
not. This makes the new method somewhat similar to the 
method of "family trees" (Bai 12003]) or longitude-time di- 
agrams (Brouwer & Zwaan [T990p . 

The literature about the longitudinal distribution of 
solar activity indicators is so wide that it is not reasonable 
to compare our results with all of them. It suffices to state 
that the general patterns revealed so far are quite similar 
to those described above. The major shortcomings of the 
previously used methods include the dependence of the 
results on some prefixed parameters or on the choice of a 
particular distribution model. 

5. Conclusions 

When introducing a new method to analyse solar activity 
patterns we started from certain methodological princi- 
ples: 

^ Input data must be homogeneous, comprehensive and 
cover as long time base as possible. 

— The analysis method must be free from any prefixed 
constants. 

— The method must not depend on the model of the ac- 
tivity indicator distributions (unimodal, bimodal etc.). 

— The computations must be as simple as possible. 

The results obtained using the new method can be ranked 
using these underlying principles. We start from the most 
evident and methodologically "clean" facts and proceed 
towards the statements which can be doubted or refined 
using additional devices. 

— The distribution of sunspots is determined by the un- 
derlying large-scale mechanism which is more persis- 
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tent than sunspots themselves. This shows up as a ten- 
dency of new sunspots to occur near the places where 
the previous sunspots were observed. 

— The mean rotation velocity of the large-scale features 
for the Northern hemisphere is 14.35 ±0.02 and for the 
Southern hemisphere 14.40 ± 0.02. 

— The rotation velocities for Northern and Southern 
hemispheres differ slightly. Consequently both veloc- 
ities manifest certain statistical averages and do not 
suggest a strong meridional coupling. 

— The large-scale patterns of activity take part in differ- 
ential rotation. The differential rotation curve is some- 
what shallower if to compare with curves obtained 
from sunspot tracking (see Zapalla & Zuccarello ll99ip . 

— The differential rotation for the Southern hemisphere 
is more similar to that obtained from sunspot analysis. 
Differential rotation curve for the Northern hemisphere 
deviates from the general rotation law especially near 
the equator. 

— The strong tendency for the spot groups to cluster 
on certain longitude dies off with time. The longest 
observable correlations can reach 15 Carrington rota- 
tions. 

— The correlations between rotations are more pro- 
nounced for lower latitudes. 

— The observation of the spot-generating mechanism be- 
ing affected by differential rotation is suggestive of 
phase mixing occurring in the solar convection zone; 
such a phenomenon is not predicted by conventional 
mean-field dynamo theory. 

The set of formulated results is impressive if to take into 
account the simplicity of the analysis method used. 
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